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ABSTRACT 

We show that the problem of constructing telescopers for 
functions of m variables is equivalent to the problem of con- 
structing telescopers for algebraic functions of m — 1 vari- 
ables and present a new algorithm to construct telescopers 
for algebraic functions of two variables. These considera- 
tions are based on analyzing the residues of the input. Ac- 
cording to experiments, the resulting algorithm for rational 
functions of three variables is faster than known algorithms, 
at least in some examples of combinatorial interest. The al- 
gorithm for algebraic functions implies a new bound on the 
order of the telescopers. 

Categories and Subject Descriptors 

1.1.2 [Computing Methodologies]: Symbolic and Alge- 
braic Manipulation — Algorithms 

General Terms 

Algorithms 

Keywords 
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1. INTRODUCTION 

The problem of creative telescoping is to find, for a given 
"function" / in several variables t\, . . . , t n , x\, . . . , x m , linear 
differential operators L involving only the ti and derivations 
with respect to the ti, and some other "functions" gi, ... , g m 
such that 

L{f) =D X1 { 91 ) + ■ ■ ■ + D Xm (g m ), 
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where D x . denotes the derivative with respect to Xj. The 
main motivation for computing such operators L (called 
"telescopers" for /) is that, under suitable technical assump- 
tions on / and the domain Q, these operators have the def- 
inite integral 




as a solution. Once differential operators for F have been 
found, other algorithms can next be used for determining 
possible closed forms, or asymptotic information, or recur- 
rence equations for the series coefficients of F. 

There are general algorithms for computing telescopers 
when the input / is holonomic [25, 15, 24, 20, 9] as well 
as special-purpose algorithms designed for restricted input 
classes [25, 26, 5]. The focus in the present paper is on two 
such restricted input classes: rational and algebraic func- 
tions of several variables. Our first result is that an algo- 
rithm for computing telescopers for rational functions of m 
variables directly leads to an algorithm for computing tele- 
scopers for algebraic functions of m — 1 variables and vice 
versa (Section 2). Our second result is a new algorithm 
for creative telescoping of algebraic functions of two vari- 
ables (Section 3), which, by the equivalence, also implies a 
new algorithm for creative telescoping of rational functions 
of three variables. The algorithm for algebraic functions is 
mainly interesting because it implies a new bound on the 
order of the telescoper in this case (Theorem 15), while the 
implied algorithm for rational functions is mainly interesting 
because at least for some examples it provides an efficient 
alternative to other methods (Section 4). 

For a precise problem description, let k be a field of char- 
acteristic zero, and k(t, x) be the field of rational functions 
in t and x = (xi, . . . , x m ) over k. Let x m denote the m — 1 
variables xi, . . . , x m -i- The algebraic closure of a field K 
will be denoted by K. The usual derivations d/dt and d/d Xi 
are denoted by D t and D Xi , respectively. Let k(t){Dt) be 
the ring of linear differential operators in t with coefficients 
in k(t). Then we are interested in the following two prob- 
lems: 

Problem 1. Given f G k{t, x), find a nonzero operator L G 
k(t){Dt) such that 

L(f) = D xi (gi) H h D Xm (g m ) for some gj G fc(t,x). 

Such an L is called a telescoper for f, and the rational func- 
tions gi, . . . , g m are called certificates of L. 



Problem 1. Given a £ k(t,ic„ 
tor L £ k(t){Dt) such that 



find a nonzero opera- 



L(a)=D xl (f}i)-\ \-D Xm _ 1 (Pm-i) for some Pj £ fc(i,x m ). 



Such an L is called a telescoper for a, and the algebraic 
functions /3i, . . . ,/3 m —l are called certificates of L. 

Both the equivalence of these two problems and the new 
algorithm for Problem 2 (when m = 2) are based on the 
general idea of eliminating residues in the input. As an 
introduction to this approach, let us consider the problem 
of finding a telescoper and certificate for a rational function 
in two variables, that is, given a rational function / £ k(t, x), 
we want to find a nonzero L £ k(t)(Dt) such that L(f) — 
D x (g) for some g £ k(t, x). We may consider / as an element 
of K(x), where K — k(t), and as such we may write 



i=l 3=1 



(1) 



where p € K [as], the Pi are the roots in K of the denominator 
of / and the otij are in K. We refer to the element on.i as 
the residue of / at /?;. Using Hermite reduction, one sees 
that a rational function h £ K(x) is of the form h = D x (g) 
for some g £ K(x) if and only if all residues of h are zero. 
Therefore to find a telescoper for / it is enough to find a 
nonzero operator L £ K{Dt) such that L(f) has only zero 
residues. For example assume that / has only simple poles, 

1. e., / = pd, b € K [x], deg x a < deg x b and b squarefree. 
We then know that the Rothstein-Trager resultant [23, 19] 

R := resultant^ (a — zD x (b), b) £ K[z] 

is a polynomial whose roots are the residues at the poles of /. 
Given a squarefree polynomial in K[z] = k(t)[z], differentia- 
tion with respect to t and elimination allow one to construct 
a nonzero linear differential operator L £ k(t)(Dt) such that 
L annihilates the roots of this polynomial. Applying L to 
each term of (1) one sees that L(f) has zero residues at each 
of its poles. Applying Hermite reduction to L(f) allows us 
to find a g such that L(f) = D x (g). 

The main idea in the method described above is that 
nonzero residues are the obstruction to being the derivative 
of a rational function and one constructs a linear operator to 
remove this obstruction. Understanding how residues form 
an obstruction to integrability and constructing linear oper- 
ators to remove this obstruction will be the guiding principal 
that motivates the results which follow. 

The authors would like to thank Barry Trager for useful 
discussions and outlining the proof of Proposition 11. 

2. TELESCOPERS FOR RATIONAL 
FUNCTIONS 

2.1 Rational and algebraic integrability 

In this section, we give a criterion which decides whether 
or not 1 is a telescoper for a rational function in fc(t,x). 
Again, let K = k(t). A rational function / £ if(x) is said to 
be rational integrable with respect to x if / = X^JLi D Xj (gj) 

for some gj £ A'(x). An algebraic function a £ K(H m ) is 
said to be algebraic integrable with re spect to x m if a = 
jyj^i D Xj -{Pi) for some /3j £ K(k m ). By taking traces, 
one can show that if a is algebraic integrable with respect 



to x m , then an antiderivative of a already exists in the 
field K(x m )(a). 

For a rational function / £ A"(x), Hermite reduction with 
respect to x m decomposes / into 



/ = A; m (r) + p 



(2) 



where r £ K{s) and a,b £ K(it m )[xm] such that deg^ (a) < 
deg Xm (b) and b is squarefree with respect to x m . It is clear 
that / is rational integrable with respect to x if and only if 
a/b in (2) is rational integrable with respect to x. Over the 
field _ftf(x m ), one can write a rational function / £ K(x.) as 

n m ; 



where p £ A(x ] and the ctij, j3i are in A"(x m ). We call 
an the Xm-residue of / at Pi, denoted by residue;c m (/, Pi). 



Proposition 3. Let f £ K (x) and P £ A"(x m ). Then 

(i) residue^ (/, P) — Q if f = D Xm (g) for some g £ A(x) 

(%%) D Xj (residue^ (/, P)) = residue^ (D x . (/), /3) for all j 
with 1 < j < m — 1. 

Proof. The first assertion follows by observing the effect 
of D Xm on each term in the partial fraction decomposition 
of g. By Hermite reduction, we can decompose / into 



Pi 

By the first assertion, either residue^ (/, /3) = a* if /9 = 
Pi or residue Xm (/, P) — if P 7^ Pi for all i — 1, ...,n. 
Applying D Xj to the two sides of the equation above yields 

D Xj {a % ) 



D x 



D *i w - E 



ctiD X] {Pi) 



E 



Then we have either residue^ (D Xj (/), P) — D Xj (cti) if p = 
Pi or residue^ (D Xj (/),/?) = if /3 ^ Pi for alii = 1, . . . , n. 
The second assertion follows. ■ 

If / is written as the form in (2), then we have 



residue^/, pi) 



D Xm (b) \ Xm =0i 



£ K{± m ){Pi). 



Therefore, all the x m -residues of / are roots of the Rothstein- 
Trager resultant (see [19, 23]) 

R := resultant^ (b, a — zD Xm (b)) £ K{k m )[z\. 

Lemma 4. Let f £ A(x). Then f is rational integrable 
with respect to x if and only if all the x m -residues of f are 
algebraic integrable with respect to x m . 

Proof. By the Hermite reduction and partial fraction de- 
composition, / can be written as 



/ = £* m w+E 



A' 



where r £ K(x), 014, Pi € K(it m ) and the Pi are pairwise 
distinct. 

Suppose that all the x m -residues on of / are algebraic 
integrable with respect to x m , i.e., on — X^jL^. 1 D Xj (ji,j) for 
some £ K(it m )(oti). Note that for each j we have 

RiihiA = Dxj (_jad_) + Dxm h^m\ . 

Then we get 

Therefore, / is rational integrable with respect to x by tak- 
ing 

n n m— 1 7-. / o \ 

^— — and g m = r + ^^ . 

1 — 1 z — 1 j — 1 

Note that all the gj and g m are in K(-x.) because ji,j £ 
K{yL m ){Pi) and Pi are roots of a polynomial in K(i. m )[xm]- 
Suppose now that / is rational integrable with respect 
to x, i.e., / = Y^JjLi D X;j (gj) for some gj £ K(x). For any i £ 
{1,2,..., n}, taking the a; m -residues of / and Y^Li D Xj (gj), 
respectively, and using Proposition 3 we get 

m — 1 

residue^/, /%) = a, = ^ D Xj (residue^ (g jt Pi)), 
3=1 

which implies that a; is algebraic integrable with respect 
to x m . ■ 

Example 5. Let f — l/(x'i +X2). Then the X2-residue of f 
at —x\ is 1. Since 1 = D X1 (xi), f is rational integrable with 
respect to Xi and X2- More precisely, 

f = D Xl f-^)+D X2 (-- ^-Y 

\Xl +X 2 J \ X±+X 2 J 

Example 6. Let f = l/(x\X2)- Then the x^-residue of f 
at is l/aii. Since l/x\ has no antiderivative in K(x\), f 
is not rational integrable with respect to x\ and x 2 . 

2.2 Equivalence 

Theorem 7. Let f £ k(t,x). Then L £ k(t)(D t ) is a 
telescoper for f if and only if L is a telescoper for every 
Xm-residue of f. 

Proof. By a similar calculation as in the proof of Proposi- 
tion 3, we have 

/.(residue^ (/, P)) = residue^ (L(f),P) (3) 

for any L £ k(t)(D t ) and p £ k(t,x m ). If L £ k(t)(D t ) is a 
telescoper for /, then L(f) = YlT=i D Xj (gj) for some gj £ 
fc(i,x). By Proposition 3 and Equation (3), for the x m - 
residue a := residue^ (/, P) at any pole P of f with respect 
to x m , we have 

m — 1 

L ( a ) = D Xj( reSidue Xm(gj,P))- 
J = l 

So L is a telescoper for a. Conversely, assume that I is a 
telescoper for any a; m -residue of /. Note that any a; m -residue 



of L(f) is of the form L(residue Im (/, /?)), which is algebraic 
integrable by assumption. Then L(f) is rational integrable 
by Lemma 4. Therefore, L is a telescoper for /. ■ 

Now we can present an explicit translation between the 
two telescoping problems by using Theorem 7. 

If we can solve Problem 2, then for a rational function / £ 
fc(i,x), first, we can perform Hermite reduction to decom- 
pose / into / = D Xm (r) + a/b; second, we compute the 
resultant R := resultant^ (a — zD Xm (b), b) £ k(t,x m )[z]; 
finally, we get a telescoper for / by constructing telescop- 
ers for all the roots of R in k(t,x m ) and taking their least 
common left multiple. 

On the other hand, if we can solve Problem 1, then for 
an algebraic function a £ fe(i,x m ) with minimal polyno- 
mial F £ k[t, x m ,a; m ], we compute a telescoper L for the 
rational function / = x m D Xm (F) / F . Note that a is the 
x m -residue of / at a. Therefore, L is a telescoper for a. 

Example 8. Consider the rational function 

2y(l - x)x(x + l)(x + 2)(t + x)(xy - y - t 4 ) 
1 - x(2 - x + (x + l)(x + 2)(t + x)(xy -y- i 4 ) 2 ) ' 

In order to find a telescoper for f , we view f as a rational 
function in y with coefficients in k(t, x) and determine its 
residues in k(t,x). Write a and b for the numerator and 
denominator of f. Since b is squarefree, the residues z of f 
are precisely the roots of the the Rothstem-Trager resultant 
resulta,nty(a—zD y (b),b) £ k(t,x)[z]. In the present example, 
these are 



x - 1 v / x(x + l)(x + 2)(x + t)' 

According to Theorem 7, it now suffices to find a telescoper 
for this algebraic function. This problem is discussed in the 
following section. 

3. TELESCOPERS FOR ALGEBRAIC 
FUNCTIONS 

We showed above how focusing on residues can yield a 
technique to find telescopers of rational functions by reduc- 
ing this question to a similar one for algebraic functions. 
In this section we describe an algorithm to solve this latter 
problem for algebraic functions of two variables. In what 
follows, the term "algebraic function" will always refer to 
functions of two variables t and x. When one tries to use 
residues to solve the problem of finding telescopers for al- 
gebraic functions one must deal with several complications. 
The first is a technical complication. One does not have a 
global way of expressing a function similar to partial frac- 
tions and so must rely on local expansions. This forces one 
to look at differentials rather than functions in order to de- 
fine the notion of residue in a manner that is independent of 
local coordinates. The second complication is a more sub- 
stantial one. There are differentials ctdx having zero residues 
everywhere that are not of the form dp — D x (P)dx, i.e. a 
is not the derivative of an algebraic function. Nonetheless, 
one knows that there must exist an operator L £ k(t)(Dt) of 
order equal to twice the genus of the curve associated to / 
such that L(a)dx — dp for some algebraic p. This will force 
us to add an additional step to find our desired telescoper. 



In Section 3.1, we will gather some facts concerning differ- 
entials in function fields of one variable that will be needed 
in our algorithm. In Section 3.2 we describe the algorithm. 

3.1 Derivations and Differentials 

In this section we review some notation and facts concern- 
ing function fields of one variable (cf. [2, 4, 8, 11, 16]). In 
the previous section the results and calculations depended 
heavily on the notion of the residue of a rational function 
of y at an algebraic function j3i of x. In the present section 
we shall also need to use the notion of a residue but since 
we are dealing with algebraic functions instead of rational 
functions, the appropriate notion is that of a residue of a 
differential uj at a place V of the associated function field E. 
We will denote this by residue-p uj and refer to the above 
mentioned books for basic definitions and properties. We 
note that when / £ E — K(x)(y), and /3i £ K(x), then 
residue y (/, /3i) = residuepaj, where uj — fdx and V is the 
place (y — j3i) of E. 

Let K be a differential field of charactersitic zero with 
derivation denoted by A (for example, K = k(t) with Dt as 
above). Let x be transcendental over K and E = K(x, y) an 
algebraic extension of K(x). We may extend the derivation 
D t to a derivation Df on K(x) by first letting Df(x) = 
and then taking the unique extension to E. We define a 
derivation D x on K(x) by letting D x be zero on K, D x (x) = 
1 and taking the unique extension of D x from K(x) to E. 
We shall also assume that the constants E Dlc = {c £ E \ 
D x (c) = 0} are precisely K. This is equivalent to saying 
that the minimal polynomial of y over K(x) is absolutely 
irreducible (cf. [10]). In [8], Chapter VI, §7, Chevalley shows 
that Df can be used to define a map (which we denote again 
by Dt) on differentials such that Df(fdx) = (Df(f))dx. 
The map Dt furthermore has the following properties: 

1. Df(dg) = d(Dfg) for any g € E, and 

2. for any place V of E and any differential uj, 

residue-p (Dfoj) = Dt (residue-p (a>)). 

Given a £ E we will want to find an operator L £ K(Df) 
and an element /3 £ E such that L(a) = D x (j3). In terms of 
differentials, this latter equation may be written as L(uj) = 
d/3, where uj = adx. 

We shall have occasion to write our field E as E = K(x, y) 
for some other x which is transcendental over K and y alge- 
braic over K(x) and work with the derivation Df defined in 
a similar manner as above. We will need to know that if we 
can find a telescoper with respect to the derivation Df then 
we can convert this into a telescoper with respect to Df . 
The following lemma and proposition allow us to do this. 

Lemma 9. Let x and x be as above and let uj be a differential 
of E . For any i = 1, 2, . . . there exists Ui £ E such that 

(ATM - (ATM = d Ui . (4) 

Proof. Write uj — adx. Lemma 1 of [16] (see also Lemma 3 
in Chapter VI, §7 of [8]) implies that 

Df{co) - Df(oj) = -d(aDf{x)). (5) 

Letting u\ = —aDf(x), we have equation (5) for i = 1. One 
can verify by induction that (5) holds for = Df(ui) — 
ViDf(x), where = (Df) 4 [aD x x] ■ D x {x). m 



Proposition 10. Let a € E, uj = adx, 

(AT + an-i(A) n_1 + • • ■ + ao 6 K{Df), 
and (3 G E such that 

((AT + an-i( Df)"' 1 + ...+ a ) (uj) = dp. 
One can effectively find f3 £ E such that 

((A)" + art-^Df)"- 1 + ... + ao )(a) = D x (/3). 
Proof. From Lemma 9 we have that 

((A)" + On-lfAfr -1 + . . . + )(W) 

= {{Df) n {uj) + du„) + a n _i((A) n_1 (^) + du n -i) 
+ . . . + a uj. 

Therefore, taking into account that the at belong to K, 

((AT + a„-i(AT _1 + • • • + ao) (w) 
= d{p — Un — a„_iu„_i - ... - aiui), 

which implies the conclusion of the proposition with f3 = 
/3 — u n — a„_iw„_i — ... — oiUi. ■ 

In the algorithm described in the next section, we will 
consider a differential uj in E = K(x, y) and assume that 

1. uj has no poles at any place above the place of K(x) 
at infinity, and 

2. the places where uj does have a pole are all unramificd 
above places of K(x). 

We describe below an algorithm that allows one to select 
an x £ E such that E = K(x,y) and such that uj satisfies 
conditions 1. and 2. above with respect to K(x). The al- 
gorithm of Section 3.2 can be used to produce a telescoper 
with respect to Df and Proposition 10 allows one to con- 
vert this telescoper to a telescoper with respect to Df . In 
the following proposition, the proof that condition 2. can be 
fulfilled was outlined to us by Barry Trager [21, 22]. 

Proposition 11. Let uj be a differential in E = K(x,y). 
One can effectively find an x £ E such that E = K(x, y) 
and 

1. uj has no poles at any place above the place of K(x) at 
infinity, and 

2. the places where uj does have a pole are all unramified 
above places of K(x). 

Proof. If 1. does not hold, let c £ K be selected so that uj 
has no poles above x = c, let 

cx 

x — . 

x — c 

This change of variables interchanges c and the point at in- 
finity, so 1. is now satisfied with respect to K(x) and we 
shall henceforth abuse notation and assume that 1. is satis- 
fied with respect to K(x). 

Let C be a nonsingular curve that is a model of E. The 
elements of E can be considered as functions on C. As noted 
in [21, p. 63], ramification occurs when the line of projection 
from the curve down to the a;-axis is tangent to the curve 
and, for each pole of uj, there are only a finite number of 
projection directions that are tangent to the curve at this 



pole. Therefore for all but finitely many choices of an inte- 
ger m, if we let x = x + my, uj will satisfy 2. with respect 
to K(x). One can refine this argument and produce a finite 
set of integers m that are to be avoided. This is done in the 
following way. 

Let M be an indeterminate and consider the field E\ — 
E(M) = ki(x,y), where ki = K{M) and x = x + My. Let 
o = if[A'f] and assume that (after a possible change of y), 
y satisfies a monic polynomial over o[x\. The behavior of 
various objects in E\ when one reduces o modulo a prime 
ideal of o is considered in [11, Chapter III, §6]. We shall 
be interested in reducing modulo ideals of the form (M — 
m), where m is an integer. One can effectively calculate an 
integral basis {wi(M)} of the integral closure of k\[x\ in E\ 
(cf. [12, 21]) and from this a complementary basis {w'^M)} 
([2, Chapter 5, §2], [4, §22]). In Chapter III §6.2 of [11], 
Eichler gives a method that will produce a finite set S C Z 
such that for m S, the set {wi(m)} is again an integral 
basis of the integral closure of K[x] in E. This method can 
be refined (and the set S slightly increased if need be) so 
that {«4(m)} is also a complementary basis. Expressing uj 
in terms of this complementary basis, 



1 



b(x 

one sees that uj will have poles precisely at the zeroes of b(x). 
If one selects m £ Z such that b(x) is relatively prime to 
D(x), the discriminant of the integral basis {wi(m)}, then 
uj will not have poles at ramification points. The finitely 
many values of m that do not satisfy this latter condition 
are roots of 

S(M) = resultants ( resultanty (b(X + MY), F(X,Y)), 
resultanty (D(X + MY), F(X, Y))), 

where F £ K[X, Y] is the minimal polynomial of y over 
K(x). m 

3.2 An Algorithm to Calculate Telescopers for 
Algebraic Functions 

We assume we are given a function field of one variable 
E = K(x, y) and a differential uj in E. We shall furthermore 
assume that uj satisfies conditions 1. and 2. of Proposition 11. 
We will describe an algorithm to find ao, ■ . . , a n £ K, not 
all zero, and P £ E such that 



(a n (D 



+ a ?l _!(AT + ■ ■ ■ + ao)(w) = dp. 



If uj = adx, then L = a n (Dt) n + a n - 1 (D't) n ' 1 + . . . + a is 
a telescoper for a with certificate p. The algorithm has two 
steps. The first step finds an operator L\ such that applying 
this operator to uj results in a differential L\(uj) with only 
zero residues. The second step finds an operator L2 of order 
at most twice the genus of E and an element (3 £ E such 
that L 2 {L 1 (uj)) = dp. 

Step 1. We will describe two methods for constructing an 
operator that annihilates the residues of uj. The first one 
requires one to calculate in algebraic extensions of K while 
the second only requires calculations in K. Throughout, let 
F(x,Y) £ _ff[a;,y] be a minimal polynomial of y over K(x) 
and let 



for some A £ K[x,y] with no finite poles and B £ K[x]. 

Method 1. We make no assumptions concerning ramification 
at the poles but for convenience we do assume that the poles 
of uj only occur at finite points. Let a £ A" be a root of B. 
For any branch of F(x, Y) — at x = a, we may write 

uj = p a (z)dz, 

where 2 = (x — a) 1 '*" for some positive integer m and p a 
is a Laurent series in 2 with coefficients in K. One can 
calculate the coefficient of I/2 in p a and this will be the 
residue of uj at this place. In this way, one can calculate 
the possible residues {ri, . . . ,r s } of uj. Let K\ be a Galois 
extension of K containing {n, . . . ,r s }. Let C be the field of 
Dt-constants in K\ and {fi, . . . ,fg\ be a C-basis of CVi + 
. . . + Cr B . Let L(Y) — wr(Y, fi, . . . re) where wr(. . .) is the 
Wronskian determinant. One sees that L(Y) is a nonzero 
linear differential polynomial with coefficients in K\ such 
that Li(rj) = for i — 1, . . . , s. Define 

Li(Y) =lclm{L CT (Y) a £ G}, 

where G is the Galois group of K\ over K, L" (Y) denotes 
the linear differential polynomial resulting from applying a 
to each coefficient of L and lclm denotes the least common 
left multiple. We then have that Li(Y) has coefficients in 
K and annihilates the residues of uj. 

Method 2. We now assume that uj has poles only at finite 
places and that there is no ramification at the poles. This 
implies that at any place corresponding to a pole, we may 
write a = X^i>i a i( x ~ x oY f° r some on £ K. Therefore 
the residue of uj at this place is 

a-i = (D x [(x 



J'd) 



This is the key to the following, parts of which in a slightly 
different form appear in [7]. 

Proposition 12. Given uj as above, one can compute a 
polynomial R £ K[Z] of degree 

m := dcg z (J?) < degy^degJB*), 

with B* the square free part of B, such that if a is a nonzero 
residue of uj then R(a) — 0. Furthermore, one can compute a 
nonzero operator L\ = a m (Dt) m +a m -i(Dt ) m_1 + . . -+ao £ 
K{Df) such that uj := L\{uj) has residue zero at all places. 



Proof. We may write 

adx = -^rdx = -^-dx + -p^rdx + 
B B± Bo 



where the A, Ai £ K(x, y) are regular at finite places and 
B = B x Bl---B\ £ K [x] is the squarefree decomposition 
of B. To achieve our goal it is therefore enough to prove the 
claim for a differential of the form adx = -^dx, where A £ 
K(x, y) is regular at finite places and B £ K\x\ is squarefree. 
Following [7], we let u be a differential indeterminate and let 



h = 



(Au 



-n\(n—l) 



(n-l)\ 



£ K{x,y)(u), 



adx 



B 



-dx 



where K(x, y){u) is the ring of differential polynomials in it 
with coefficients in K(x, y) and (. ..)'*' denotes i-fold differ- 
entiation with respect to x. Let V be a place where a has 
a pole and let a and b denote the values of x and y at the 
place. We note that since A is regular at V and V is not 



ramified, any derivative of A is also regular at V (one needs 
the hypothesis that these places are unramified to make this 
claim). Taking into account the rules of differentiation, we 
see that 



h 



(x,y,u,u 



q(x)u t 



where p(x, Y, z , zi, ... , z„-i) £ K[x, Y, z , Zi, ■ . . , z n -i], t is 
some positive integer and q(x) £ K[x\ does not vanish at V, 
i.e. q(a) 7^ 0. Let 



and 



V = p(x, Y, B , iS", ±B^,. . . , ±B W ) £ K[x, Y] 



q = q (x)(B'Y eK[x] 



One then shows, as in [7], that p(a, b)/q(a) is the residue of 
■j^dx at V. 

The above argument shows that the polynomial 

R = resultant^ (resultant y (p - Zq, F),B) £ K[Z] 

vanishes at the residues of adx. The degree estimate for R 
follows from the general degree estimate for resultants which 
states for any S,T £ K[u, v] that deg u (resultant^ (5, T)) is 
at most 

degJS) deg„(T) + de gv (S) deg„(T). 

This implies first that the inner resultant in the definition 
of R has Z-degree at most deg Y (F). (Note that no degree 
estimates for p and q are needed because deg z (F) — 0.) Ap- 
plying the rule again to the outer resultant gives the desired 
bound deg Y (F) deg x (B). 

Let R £ K[Z] be the polynomial above. If necessary, 
we may replace R by a squarefree polynomial having the 
same nonzero roots so we shall assume that R is square- 
free and of degree m. Using the fact that R and 45 are 
relatively prime, there exist polynomials Ri £ K[Z] of de- 
gree at most m — 1 such that if 7 is a root of R, then 
Dl(y) = Riiy) for i = 0, 1, . . .. Since each Ri has degree 
at most m — 1, there exist a m , . . . , ao £ K, not all zero, 
such that (a m (D x t ) m + am-i(i)i)'"^ + . . . + a () )( 7 ) = 
for any root 7 of R. Using the fact that residue-p (DtLj) = 
Df (residue-p(tj)) for any place V, one sees that for L\ = 
a m (D?) m + a m -i(Df) m_1 + ■ • • + ao, w = L x (oj) has zero 
residue at any place. ■ 

Remark. Although Method 2 does not require calculations in 
an algebraic extension of K , one needs the condition on ram- 
ification to prove that it is correct. This condition is painful 
to verify and although Propositions 10 and 11 imply that we 
can make a transformation, if necessary, to guarantee that 
the differential has poles at places that are not ramified, 
making such a transformation can increase the complexity 
of the data. In practice, one should calculate the operator L\ 
above without testing if the places at poles are ramified, cal- 
culate the operator L2 as in step 2 below (which requires no 
assumption concerning ramification) and then test to see if 
the resulting operator L 2 o L\ is a telescoper by checking if 
the identity Li(L\{pi)) = D x (f3) holds, a simple calculation 
in K(x, y). If this equality does not hold, then one can make 
a change of variable x := x + my for a random m and try 
again. Proposition 11 guarantees that after a finite number 
of trials one will succeed. 



Example 13 (continuing Ex. 8). Let F = y 2 — x(x + l)(x + 
2)(x + t) and consider 

( t 4 1 \ , u , 
uj — 1 — )dx = —ax. 

\ x — 1 y) v 

where u — (x — l)y + t 4 x(x + l)(x + 2)(t + x) and v — 
x(x + l)(x + 2)(x + t)(x — 1) . The only pole of ui is a simple 
pole at x — 1, so the residues of u) are the roots of 

residue x (residue y (u — zD x (v), F),v) = (. . .)(z — t 4 ) 2 z 8 , 

where (. . .) stands for some factors which are free of z and 
therefore irrelevant here. The only nonzero residue t 4 is an- 
nihilated by Li := tDt — 4, so 



Cj = (tD t -4)(w) = 



has no nonzero residues. 



(9t + 8x)y 



2x(x+l)(x + 2)(x + t) 2 



dx 



Remark. 1. In [21], Trager develops a Hermite reduction 
method for algebraic functions which, when applied 
to the differential uj above, shows how one can write 
w = (Di(ji) + g2)dx, where g\,g% £ E and (72 has 
only simple poles at finite points. Regretably, g% may 
have poles (of higher order) at infinity. Nonetheless, it 
would be interesting to see if Trager's procedure can 
be used to increase efficiency in our algorithm. 

2. The above argument strongly relies on the fact that 
we are assuming that the places where ui has poles 
are not ramified above places in K(x). It would be 
of interest to give a method to calculate an operator 
L\ satisfying the conclusion of Proposition 12 without 
this assumption. 

Step 2. Let ui be as in the conclusion of Proposition 12. 
Again using the fact that residue-p(Dfd>) = Df (residue-p(iZ))) 
for any place V, we have for all i £ Z that (Df) l (uj) is again 
a differential with zero residues at all places. Such a differ- 
ential is called a differential of the second kind ([8], p. 50) 
and a differential of the form dy, 7 £ E is called an exact 
differential. Note that any exact differential is a differen- 
tial of the second kind. Corollary 1 of ([8], p. 130) states 
that the factor space of the space of differentials of the first 
kind by the space of exact differentials is a if-vector space 
of dimension equal to 2G, where G is the genus of E. There- 
fore, there exist 0,20, ■ ■ ■ ,0,0 £ K, not all zero, such that for 
L 2 = &2g(D? ) 2G + ^c-itA 1 ) 20 " 1 + .. . + So, L 2 (Cj) = d(3 
for some fJ £ E. Such Li and f3 can be found as follows. 

Let u — adx and let [E : K(x)\ — m. For each i > 0, 
there exist ai,o, ■ • ■ , cti.m-i £ K(x) such that 



(Dty(a) = (y,...,y m - 1 ) 



( «i,0 



In addition, there exists an m x m matrix A with entries in 
K(x) such that 

(D x {y),...,D x (y m - 1 )) = (y,...,y m - 1 )A. 

Let ao, . . . , 02G be elements of K and f3o, . . . , f3 m -i elements 
of K(x). Letting f3 = /3 +/3ij/+. . .+(3 m -iy m ~ 1 , the equation 

d/3 = (a2G(Dt) 2G '(&) + . . . + aoa)dx 



is equivalent to 
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In [3] , Barkatou describes a decision procedure for deciding if 
there exist nontrivial j3o, . . . , Pm-i & K(x) and ao, . . . , a^G £ 
K satisfying (6) when K is a computable field (i.e., the arith- 
metic operations and derivation are computable and one has 
an algorithm to factor polynomials over K). Therefore one 
can apply this to K = k(t), where k is a computable field of 
characteristic zero to produce a desired L2 and /3. 

Example 14 (continuing Ex. 13). Let again F = y 2 — x(x + 
l)(x + 2)(x + t) and consider the differential 

(9t + 8x)y 

UJ — ; — — rrrdx. 

2x(x + l)(x + 2)(x + t) 2 

Since the field E has genus 1 and u> has only zero residues, 
there exists a telescoper for Q of order 2. Indeed, the algo- 
rithm outlined above finds that L2{ui) = dp, where 

L 2 = 4(99t 5 - 540t 4 + 1055* 3 - 870t 2 + 256t)D 2 

+ 4(297t 4 - 1269t 3 + 1900t 2 - 1152* + 256)A 

+ 3(99t 3 - 306t 2 + 307t - 96) 

and 

_ 3(429t 3 + 330i 2 a; - 89H 2 - 648ta: + 384t + 256x)y 
' (t + xj 3 ' 

For the differential ui from Example 13, it follows that we 
have Luj = d/3 with 

L = L 2 o (tD t - 4) = 4(t - 2)(i - l)t 2 (99t 2 - 243* + 128) A 3 

+ U{99t 4 - 189i 3 - 210i 2 + 588t - 256)D 2 

- 3(1089t 4 - 4770t 3 + 7293t 2 - 4512f + 1024)D t 

- 12(99t 3 - 306f 2 + 307t - 96). 

By Theorem 7, this operator L is also a telescoper for the 
trivariate rational function f from Example 8. Certificates 
g, h with 

L(f) = D x (g) + D v (h) 

can be obtained from (3 following the calculations in the proof 
of Lemma 4- They are however too long to be printed here. 

Remark. Telescopers and certificates for holomorphic differ- 
entials arise in Manin's solution of Mordell's Conjecture [16, 
17] and Step 2 of our procedure is just an effective version 
of considerations that appear in these papers. Telescopers 
for holomorphic differentials are also referred to as Gauss- 
Manin Connections. 

Combining the estimates on the order of the operators 
computed in steps 1 and 2 gives the following bound on the 
order of telescopers for algebraic functions. It can be viewed 
as a generalization of Corollary 14 in [5], which says that 
for every rational function / = A/B £ K(x) there exists a 
telescoper of order at most deg^. B* , where B* is the square 
free part of B. 



Theorem 15. Let E be an algebraic extension of K(x), a — 
A I B 6 E so that A is regular at finite places and B £ K [x]. 
Let B* be the square free part of B. Then there exists /? £ E 
and a nonzero operator L £ K{Dt) with L(a) = D x (fi) and 

deg Dt (i) < [E : K(x)]deg x (B*) + 2genus(E). 

4. IMPLEMENTATION AND OTHER 
EXAMPLES 

We have produced a prototype implementation of the al- 
gorithms described above on top of Koutschan's Mathemat- 
ica package "HolonomicFunctions.nl" [13] and compared the 
performance to the built-in creative telescoping implemen- 
tations of this package. In order to make the comparison as 
fair as possible, we have tried to reuse as much code from 
Koutschan's package as possible, so that the timings will not 
implicitly compare two different implementations of some 
subroutine but reflect as closely as possible the speed-up (or 
slow-down) offered by the ideas presented above. 

Five different methods to solve the creative telescoping 
problem for a rational function / £ k(t,x,y) were consid- 
ered: (CC) first use Chyzak's algorithm [9] to find a holo- 
nomic system S of operators in k(t, x){Dt, D x ) such that 
for all L 6 S there exists a rational function g £ k(t, x, y) 
with L(f) = D y (g), afterwards apply the same algorithm 
to S to obtain a telescoper L £ k(t)(Dt) for /; (CK) first 
compute S C k(t, x)(Dt, D x ) as in variant (CC), then apply 
Koutschan's ansatz [14] to 5* to obtain a telescoper L for /; 
(K) compute a telescoper for / directly with Koutschan's 
ansatz; (EC) use the reduction from Section 2, then ap- 
ply Chyzak's algorithm to the resulting algebraic functions, 
and then take the least common left multiple of the results; 
(EA) use the reduction from Section 2, then apply the al- 
gorithm from Section 3 to the resulting algebraic functions, 
and then take the least common left multiple of the results. 

Table 1 shows the performance of these five approaches 
for the following examples. 

1. The rational function / from Example 8 above. This 
example is not representative but was particularly de- 
signed to be easy for our algorithms and difficult for 
the known ones. 

2. Here / := ±h{±,x,y) with h(t,x,y) = (l - ^ - 

1 — y 1— t 1 — xy 1 — xt 1 — yt 1 — xyt/ 

the problem of enumerating diagonal 3D-Queens walks 
raised in [6]. Our calculation confirms the correctness 
of the telescoper conjectured there. 

3. Let now h(t,x,y) = (l--^ - -^- t - j^-) _1 

v ' '^z V 1 — xy 1 — xt 1—yt 1 — xyt/ 

and f — —h(-k-,x,y). This is a variation of the 

J xy V x y ' ' " ' 

previous problem, with the points (2n, n, n) replacing 
the diagonal and now allowing steps along the axes. 

4. The rational function h(t,x,y) = 2t 2 /((l - i)(3 - (x + 
y+t + xy + xt + yt) + 3xyt)) appears in [18] as the gener- 
ating function for the probability of certain structures 
in random groves. See [18] for details on the combi- 
natorial background. Here we compute the diagonal 
series coefficients of / by applying creative telescoping 
t° / = xy^i^^^ x i y)- ^ s can ^ e seen in this example, 
our algorithms are in some cases not superior. 
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CK 
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EC 


EA 
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order 


degree 


bytecount 


1 

1 


s> 1UUI1 




luy.uo 


l.OU 
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3464 
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16029.55 


40043.01 


>100h 


1390.14 


1646.53 
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71 


76472 


3 


>150h 


350495.88 


>150h 


203.44 


328.08 


9 


93 


140520 


4 


638.70 


1099.08 


>40Gb 


37606.28 


216201.88 


10 


32 


41840 


5 


23823.70 


676.13 


19085.67 


1114.34 


3117.43 


7 


27 


25320 



Table 1: Runtime comparison for the examples described in the text. 



5. With h as before, we now consider / 



-h(- 



,x,y). 



[12 



Note the large difference between CC and CK. 



We have put timings for a number of additional examples [13 
on the website [1]. Also our code and the certificates for 
Example 13 can be found there. The examples we tested 
suggest that the reduction from rational functions to alge- [14 
braic functions can cause a decent speed-up. It does seem 
to depend on whether the Rothstein-Trager resultant of the 
input factors into several small factors or not. If it does, [15 
it is advantageous because solving several small instances of 
Problem 2 is cheaper than solving a single big one. Whether [16 
after the reduction, the algorithm of Section 3 or some other 
method is applied to the resulting algebraic functions, makes 
usually not much of a difference. Our algorithm tends to be 
faster when Step 1 in Section 3.2 already finds a great part [17 
of the telescoper, leaving only a small coupled differential 
system to be solved in Step 2. 
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